Asymptotic behavior of an elastic beam fixed on a small 
part of one of its extremities 

(Preliminary version of a Note to be published in a 
slightly abbreviated form in 

C. R. Acad. Sci. Paris, Ser. 1, 338 (2004), pp. 975-980) 

Juan CASADO-DIAZ 3 , Manuel LUNA-LAYNEZ 3 , Francois MURAT b 

a Departamento de Ecuaciones Diferenciales y Analisis Numerico, Universidad de Sevilla, c/ Tarfla s/n, 41012 
Sevilla, Espana 

b Laboratoire Jacques-Louis Lions, Universite Pierre et Marie Curie, boite courrier 187, 75252 Paris cedex 05, 
France 



Abstract. We study the asymptotic behavior of the solution of an anisotropic, heterogeneous, linea- 
rized elasticity problem in a cylinder whose diameter e tends to zero. The cylinder is assu- 
med to be fixed (homogeneous Dirichlet boundary condition) on the whole of one of its 
extremities, but only on a small part (of size er E ) of the second one; the Neumann boundary 
condition is assumed on the remainder of the boundary. We show that the result depends 
on r £ , and that there are 3 critical sizes, namely r E = e 3 , r e = e, and r E = e 1 ^ 3 , and in 
total 7 different regimes. We also prove a corrector result for each behavior of r e . © XXXX 
Academie des sciences/Editions scientifiques et medicales Elsevier SAS 

Comportement asymptotique d'une poutre elastique fixee sur une petite partie 
de Vune de ses extremites 

Resume. Nous etudions le comportement asymptotique de la solution d'un probleme d'elasticite 
lineaire anisotrope et heterogene dans un cylindre dont le diametre e tend vers zero. Le 
cylindre est fixe ( condition de Dirichlet homogene ) sur la totalite de I 'une de ses extremites, 
mais seulement sur une petite partie {de taille er E ) de V autre base; sur le reste de la 
frontiere on a la condition de Neumann. Nous montrons que le resultat depend de r e , et 
qu'il existe 3 tailles critiques, a savoir r £ = e 3 , r e = e et r e = e 1 ^ 3 , et au total 7 com- 
portements differents. Nous donnons un resultat de correcteur pour tons les comportements 
de r e . © XXXX Academie des sciences/Editions scientifiques et medicales Elsevier SAS 



Version francaise abregee 

Dans cette Note, nous etudions le comportement asymptotique de la solution d'un probleme d'elasticite 
lineaire anisotrope et heterogene pose dans un cylindre fl £ — (0, 1)xe5c K 3 dont le diametre e tend vers 
zero et dont l'axe est le premier axe de coordonnees Ox\. A l'une de ses extremites {x\ — 1) le cylindre 
est fixe (condition de Dirichlet homogene) sur la totalite de la base Y\ = {1} x eS, tandis qu'a 1' autre 
extremite (x\ — 0), il est seulement fixe sur une petite partie Tq = {0} x er e So de la base, partie qui 
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est bien plus petite que e car r £ tend vers zero. Sur le reste de la frontiere de fl £ , on a la condition de 
Neumann (voir Figure 1). Pour A tenseur d'ordre 4 coercif a coefficients continus dans VL{VL = (0, 1) x S), 
et f e L 2 (Vif, h G L 2 (fl) 3x3 donnes, on definit A £ , F £ et H £ par Q et et on definit U £ comme la 
solution du probleme d'elasticite Q. 

Le but de cette Note est de decrire le comportement asymptotique de U £ et de donner un resultat de 
correcteur pour e(U £ ) quand e et r £ tendent vers zero. Ce resultat est decrit dans le Theoreme enonce dans 
la version anglaise ci-dessous, qui fait apparaitre 3 tailles critiques, a savoir r £ w e 3 , r £ » e et r £ w e 1 / 3 , 
qui separent 4 zones correspondant ar E < e 3 , e 3 « r £ C £, e « r £ « e 1 / 3 et e 1 / 3 <C r £ ^ C, done 
au total 7 cas differents. Le resultat de convergence est donne par oil it est defini par la solution de ©, 
tandis que (jfji est un resultat de correcteur pour e(U £ ). 

Quand r £ <C e 3 , l'ensemble Tq = {0} x er £ So est si petit que la condition de Dirichlet imposee a 
U £ pour x\ = disparait completement a la limite. Quand e 3 <C r £ <C £, l'ensemble Tq est suffisamment 
grand pour qu'ala limite on ait Ca(0) = Opoura G {2,3}.Quande <C r £ <C e 1 / 3 , on a en outre Ci(0) = 0. 
Finalement, quand e 1 / 3 -C r £ , l'ensemble Tq est si grand que toutes les conditions de Dirichlet possibles 
sont satisfaites a la limite, e'est a dire Ca(0) = Ci(0) = f^(0) = c(0) = (voir la version anglaise pour 
les definitions des espaces fonctionnels qui interviennent dans le probleme limite @ et pour les definitions 
de ces fonctions). 

Sauf quand r £ est de taille critique (e'est a dire quand r £ w e A avec A = 3, 1 ou 1/3), la forme bilineaire 
B qui intervient dans © et la fonction P £ qui intervient dans l|6) sont nulles. Mais pour ces trois tailles 
critiques, B est une forme bilineaire coercive (en Cq(0) pour A = 3, en Ci(0) pour A = 1, et en c(0) et 
(0) pour A = 1/3, voir Q de la version anglaise), de sorte que pour ces tailles critiques r £ w e A , B est 

une penalisation, par rapport au cas r £ <C e x , des nouvelles conditions de Dirichlet qui apparaitront pour 
r £ 3> £ A , ce qui retablit une certaine continuite dans la transition. 

La presente Note est la generalisation au cas de l'elasticite lineaire de 1' etude menee dans |2|, |3| pour 
le probleme de diffusion analogue quand le parametre t £ qui intervient dans ces travaux est t £ = 0. Les 
demonstrations detaillees seront donnees dans |4|. 

1. Position of the problem and notation 

In this Note we study the asymptotic behavior of the solution of an anisotropic, heterogeneous, linearized 
elasticity problem posed in a thin cylinder fl £ whose diameter e tends to zero and whose axis is the first axis 
of coordinates Ox\. On one of its extremities (xi = 1) the cylinder is fixed on its whole basis Tf whereas 
on the second one (xi — 0) it is fixed only on a small part Tq of it, of diameter er £ much smaller than e. 
The Neumann boundary condition is assumed on the remainder of the boundary of il £ . Mathematically the 
problem can be formulated as follows. 

For e > 0, we consider r £ a positive parameter which tends to zero with e. Let and S be two bounded 
smooth domains in M 2 , with G S. We define fi = (0, 1) xS,Q £ = (0, 1) x eS and r £ = T £ U Tf , where 
Tq = {0} x er £ So, T\ = {1} x eS. Observe that the size of Tq is much smaller than the size of the basis 
{0} x eS since r £ tends to zero. The thin cylinder Q, £ is represented in Figure 1 in the case where So and 
S are both balls of M 2 . 
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An elastic beam fixed on a small part of one of its extremities 

The elements of R 3 are decomposed as x = {x\,x'), x\ £ R, x' = (2:2,2:3) 6 R 2 . We denote by 
{e 1 , e 2 , e 3 } the canonical basis of R 3 and by £(R 3x3 ) the space of linear maps of R 3x3 into itself (or in 
other terms of fourth order tensors), where R 3x3 is the space of the 3x3 symmetric matrices. We adopt 
Einstein's convention of repeated indices. Greek indices (a and (J) take only the values 2 and 3, while latin 
indices (i and j) take the values 1, 2 and 3. 

We consider A 6 C°(Tl; £(R 3x3 )) such that there exists m > with 

A(y)ti^m\t\ 2 , V£eR 3x3 , VyeTl, 
and we define A £ e C°(f7 £ ; £(R 3x3 )) by 

A £ (x) = A(xi,—), Vxetf. (1) 

e 

We also consider / £ L 2 {Uf and h e L 2 (f2) 3x3 , and we define F £ € L 2 {U £ ) 3 and H £ e L 2 (Q £ ) 3x3 by 
F e (x) = f 1 (x u -)e 1 + ef a (x 1 ,-)e a , H £ (x) = h(x u -), a.e. x e £ . (2) 

£ £ £ 

In the thin domain fi e we consider the elasticity problem 

r U £ e H^(n £ ) 3 , 

I f A e e(U e ):e(U e )dx= f F'TI'dx+f H E :e{U £ )dx, VU'eH^in 6 ) 3 , <J) 

where 

H^ E (n e ) = {U e H 1 ^) : U = 0onT £ }; 

observe that in the above formulation, as well as in the remainder of the present Note, complex numbers 
never appear, and that U (and later u, v, W) denotes the test function associated to the solution U £ , and 
not its complex conjugate. Observe also that the solution U e of Q satisfies a non homogeneous Neumann 
boundary condition on the part dfl £ \ T £ where it is not fixed, since integrating by parts J nE H £ : e(U £ )dx 
(when h and therefore H £ is sufficiently smooth) produces both body forces and surface forces. Similarly to 
the body forces F £ we could have introduced explicit surface forces G £ on d$l £ \ T £ , but we have preferred 
not to include them for the sake of simplicity. 

It is well known that problem Q has an unique solution (see, e.g., |5|). The aim of the present Note, 
which announces our paper |4|, is to describe the asymptotic behavior of the solution U £ and to give a 
corrector result for e(U £ ) as e tends to zero. The result depends on the behavior of r £ and exhibits 3 critical 
sizes, namely e 3 , e, and e 1 / 3 , so that there are 7 different regimes: r £ -C e 3 , r £ e 3 , £ 3 « r E C e, r £ w e, 
£ <C r £ <C e 1 / 3 , r £ e 1 / 3 , and e 1 / 3 « r e ^ C, where r £ <C £ X stands for r e /e A — -> (and equivalently 
£ X <C r £ for r £ /e A — * +00), while r £ « e A stands for r £ /e A — > p, for some p with < p < +00. 

To express the results and to make the proofs, we will use two changes of variables. The first change of 
variables is given by 

x' 

y = y £ {x) with y 1 =x 1 ,y' = —, 

£ 

which transforms the variable domain il £ into the fixed domain ft. This is the usual change of variables 
used to study equations in thin cylinders (see, e.g., |9|, 1 10|, 1 11 1). When r £ — 1 and Sq — S (but 
the same proof works for r £ = C independent of e such that CSq C S), it was used successfully in |6|, 
l9l . [10| to pass to the limit in (|3}- But when r £ tends to zero with e, this first change of variables does not 
provide the information we need about the behavior of U £ in the part of £1 £ close to Tq. Thus we introduce 
a second change of variables given by 

z — z £ (x) with z= , 

£T £ 

which transforms the variable domain il £ into a variable domain Z £ , the limit of which is the half space Z = 
(0, +00) x R 2 . Observe that the Dirichlet boundary condition is now imposed on the fixed part {0} x Sq 



3 



J. Casado-Diaz et al. 



of the boundary of Z e . This change of variables provides a suitable rescaling near X\ — 0. It was used 
successfully in |2|, to study the diffusion problem similar to in the geometry that we consider 
in the present Note (and even in a more complicated one, where ft 5 is made of union of two cylinders 
{(— < e ,0) x er £ S'o}U{(0, 1) x eS} and where the Dirichlet condition is imposed on x\ = — i e andxi = 1; 
the geometry considered in the present Note corresponds to t £ = 0; a problem of conduction in a notched 
beam of the same type was solved in 1 1 1 by using the same change of variables). 
We denote by D 1,2 (Z) the Deny's space 

D 1 > 2 (Z) = {p:peL*(Z), VpeL 2 (Z) 3 }. 

We will also use the functional spaces (already used in 1 7 1 , 1 8 1 , |9|, 1 1 1 ) 



BN, 



"&(«) = {«: 3C« Gi/ 2 (0,1), Ca(l) = ^(l) = 0, u a (y) = Ca(yi), Va£{2,3}, 



3d e H\o, l), Ci(i) = 0, = Ci(y) - ^(vi)v«}, 

R b (fl) = \v: VxGL 2 (0,l;H\S)), [ Vx( yi , y')dy' = a.e. Vl G (0, 1), 
< 1 Js 

3c e H\0, 1), c(l) - 0, v 2 (y) = c( yi )y 3 , v 3 (y) = -c( yi )y 2 }, 
RD 2 L {Q) = {w: m =0, w a € L 2 (0,1; H\S)), f w a (y u y')dy' = 0, Va G {2, 3}, 
{V3W 2 {yi,y') -y 2 w 3 (y 1 ,y'))dy'= a.e. t/i G (0,1) j. 



For a given (u, v, w) G BNb(Si) x i?&(f2) x (17), we denote by -E(u, u, io) the second order symmetric 
tensor with values in K;? x3 defined by 

E n (u,v,w) = e n (u), Eip(u,v,w) = eip(v), E a p(u,v,w) = e a p(w), Va,j3 G {2,3}. 

2. The result and some comments 

The asymptotic behavior of the solution of Q depends on the size of r e with respect to e. Seven regimes 
appear in the following Theorem which describes the asymptotic behavior of II s and provides a corrector 
result for U e and e(U e ). 

THEOREM — Let U £ be the solution o/@. There exist a closed linear subspace £ of BN}, (Q) X Ri, (Q) x 
RD^i^l), a function P e G L 2 (fl £ ) 3x3 , and a bilinear continuous nonnegative form B on E x £ such that, 
defining (it, v, w) as the solution of the variational problem 

(u, v, w) G £, 

AE(u, v, w) : E(u, v, w)dy + { B(u, v, w), (u, v,w)) = 



(4) 



n 



then, when e tends to zero, we have 



fudy+ / h : E(u,v,w)dy 7 V(u,v,w) £ £, 



ttL / ( \U{{x) - Ul (x u ^)| 2 + E \ eU &) 

' ' • /Q ' V Q=2 



u a (xx)\ 



dx 



1 



\e(U e )(x) - E(u,v,w)(x u 



pc( — )\idx 



(5) 
(6) 
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The definitions of £, P £ and B do not depend on the forces f and h which define F £ and H £ , but only on 
the set So, on the fourth order tensor A, and on the behavior ofr e when e tends to zero, and more specifically 
of its behavior with respect to £ 3 , e, and e 1 / 3 , such that there are 7 regimes, which are described now. 

• Ifr e < e 3 , then £ = BN h (Q) x R b (Q) x RD£(Q), P £ = 0, B = 0. 

• Ifr £ « e 3 with r £ /e 3 -> p, then £ = BN b (Vl) x R b (fl) x RD£(Cl), and defining <p a , a G {2,3}, 
as the solution of 

tp a G D 1,2 (Z) 3 , <p a = e a on {0} x S , 

[ A(0)e(tp a ) : e{lp)dz = 0, Vtp E D 1 ' 2 {Zf, Tp = on {0} x S , 

then one has 

P £ {z) = -^Ca(0) e(<p a )(z), a.e. z E Z, 
(B(u,v,w),(u,v,w))=pj z A(0) (Co (0) e(<p a )) : (0) e(</)) dz, V(u, «,H(«,P)e«. 

• /fe 3 <C r £ <C e, ?nen 

£ = {(u,«,u;) G SJV 6 (fi) x R b (fl) x RD^(n) : ( a (Q) = 0, Va G {2,3}}, P £ = 0, 6 = 0. 

• Ifr £ w £ vvz'fn r £ /e — > p, f/zen 
£ = {(u,v,w) E SiV 6 (f7) x R b (Q) x RD^(n) : ( a (0) = 0, Va G {2,3}} , 

ant/ defining if 1 as the solution of 

tp 1 G D 12 (Zf, tp 1 = e 1 on {0} x S , 

J ^(O)e(^ 1 ) : e(lp)dz = 0, \/p E D 12 {Z f, Tp = on {0} x 5 , 

and then setting ip 1 = (p 1 + a a tp a , where (a 2 , a^) is defined by 
(a 2 ,a 3 ) GR 2 , 

y A(0)(e(^) + a«e(^)):(«/9e(/))^ = 0, V(a 2 ,a 3 ) G M 2 , 

fnen one was 

P £ (z) = --^Ci(0)e(^)(z), a.e. z E Z, 
{B{u,v,w),(u,v,w))=pJa(0) (Ci(O)e(^ 1 )) : (Ci (0) e(<p 1 )) tfz, V(u, u, w), (u,v,w) G £. 
//e <C r £ <C e 1/3 , f«en 

£ = {(u,w,w)GB7V b (fi)xi?,,(fi)xi?L> 2 L (fi) : ^ a (0)=Ci(0)=0, VaG{2,3}} , P £ = 0, B = 0. 

//> £ w e 1 / 3 vw'rn (r £ ) 3 /e — ► p, f/zen 

£ = {(u,v,w) G BA^ft) x i? b (fi) x RD^(Q) : ( a (0) = Ci(0) = 0, Va G {2,3}} , 
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and defining if) 1 as the solution of 

i) 1 G D 1 ' 2 (Z) 3 , tP 1 = z 3 e 2 - z 2 e 3 on {0} x So, 

A{Q)e{^) : e($)dz = 0, G D^ 2 {Z) 3 , ip = on {0} x S , 



one/ i/) Q , a G {2, 3}, as the solution of 

t[, a G D 1,2 (Z) 3 , t/> q = zie" - zae 1 on {0} x S , 

A{ti)e{ip a ) : e(^)dz = 0, Vip G D 1,2 (Z) 3 , ^ = on {0} x S , 



anaf fnen setting Tp 1 — ip l + 6^, where (b\, b 2 , b\), i G {1, 2, 3}, is defined by 
(6i,6|,4)GK 3 , 

A(0) (e(^) + 6le(/)) : (Sje(^)) dz = 0, V&XX) 



then one has 



P £ (z) = ~ (c(0) e(^)(z) + ^(0) e$ a )(z)) , a.e. z G Z, 



(B(u,v,w),(u,v,w))=p [ A(0) f C (0) e (^) + ^(0)e(^; 

: ^c(0)e(^ 1 ) + ^(0)e(^ a )j dz, V(u, v,w), (u, v, w) G £. 
Ife 1/3 < r e < C, then 

£ = Uu : v,w)eBN b (Q)xR b (n)xRD 2 L (n) : ( a (0)= fi(0)= ^L(0)= c(Q)= 0, Vae{2,3}|, 
I dyi J 



P £ = 0, 6 = 0. 



Let us make some comments about the statement of this Theorem. 

Assertion (jfji is a corrector result, since using U = U e as test function in and Korn's inequality, one 
can prove that 



1 



\e{U £ )\ 2 dx < C 



1 

W\ 



\f\ 2 dx- 



1 



\h\ 2 dx 



and since one can also prove that 
1 



x\, 2 I 



\E(u,v,w)(xi, — )| da; + 



|P £ ( 1. 

1 w 



If one examines the definition of £, one realizes that the number of Dirichlet boundary conditions imposed 
in the definition of £ increases with the size of r £ . Indeed, in view of the definitions of BNb(£l), Rb(ty 
and RD 2 (fl), the sole functions which have a trace for xi = are £ a , £i, and c, where a G {2,3} (the 
other functions, namely v\ and u> a , have no trace for x\ = 0). When r e <C £ 3 , the set Tq = {0} x er e So is 
too small and the homogeneous Dirichlet boundary condition imposed for x\ = to the solution U e of Q 
completely disappears at the limit. When e 3 -C r e <C £, the set Tq is sufficiently large to impose at the 
limit that £ Q (0) = for a £ {2,3}. When e < r £ < £ 1/3 , one further has Ci(0) = 0. Finally when 
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An elastic beam fixed on a small part of one of its extremities 

£ i/3 ^ r e ^ ^ set pe j s so j ar g e that all the possible Dirichlet boundary conditions are imposed at x\ = 0, 
namely G (0) = Ci (0) = §f (0) = c(0) - 0. 

Except in the three regimes where the size of r £ is critical (i.e. when r £ w £ A , with A = 3, 1, or 1/3), 
one always has P £ — and B = 0. For these three critical sizes, one can show that B is a coercive bilinear 
form, in the sense that there exists some n > such that 

(B(u,v,w),(u,v,w)) ^ np \_. ICa(0)| 2 , whenr 6 «£ 3 , 

a=2,3 

i (B(u,v,w),(u,v,w)) ^ np|Ci(0)| 2 , when r £ w e, ^ 

(£(«,«,«;),(«, >np( |c(0)| 2 + £ 

V a=2,3 

This implies that for every critical size r £ ss e A , the new Dirichlet boundary conditions which appear for 
r e 3> £ A (with respect to those imposed for r £ <C £ A ) are penalized by the value of p. This introduces some 
type of continuity in the transition of the Dirichlet condition between the two regimes which are separated 
by a critical size e A . For these three critical sizes, the functions ip a , ip 1 , ip 1 , %p l , and ip % are in some sense 
generalized capacitary potentials of {0} x So in Z, and the bilinear form B is in some sense an asymptotic 
trace of some type of capacity of Tq in fl £ for the weighted energy J Qe A(x)e(<p) : e(tp)dx. 

The present work is the natural generalization to the elastic case of 0, 0, where diffusion problems 
were posed in the union of two cylinders {(— 1 £ , 0) x er e S'o}U{(0, 1) x eS}, with both t £ and r £ tending to 
zero (the present geometry corresponds to t £ = 0). When t £ = 0, the diffusion problem was in comparison 
more simple, since only one critical size, namely r £ « e, appeared in the analysis, separating a pure 
Neumann boundary condition [corresponding to the analogue of u satisfying u G i? 1 (0, 1), u(l) = 0] 
for r £ <C £ and a pure Dirichlet boundary condition [corresponding to the analogue of u satisfying u G 
i? 1 (0, 1), u(0) = u(l) = 0] for r £ » e. These works were related to Cl> where a notched beam for 
diffusion problems was considered. The present work is also related to 1 7 1, 1 8 1, where a multidomain made 
of an elastic vertical beam of length 1 and of radius r £ and of an horizontal plate of radius 1 and of height e 
was considered. 

The detailed proofs of the results of the present Note will be given in a forthcoming paper |4|. 
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when r £ 
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